Electron Waiting Times in Non-Markovian Quantum Transport 
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We formulate a quantum theory of electron waiting time distributions for charge transport in nano- 
structures described by non-Markovian generalized master equations. We illustrate our method by 
calculating the waiting time distribution of electron transport through a dissipative double quantum 
dot, where memory effects are present due to a strongly coupled heat bath. We consider the influence 
of non-Markovian dephasing on the distribution of electron waiting times and discuss how spectral 
properties of the heat bath may be detected through measurements of the electron waiting time. 
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Introduction. — The theory of open quantum systems 
is important in many branches of physics [1] . The typical 
scenario consists of a small quantum system with a few 
degrees of freedom coupled to a large environment. The 
system evolves coherently due to its internal dynamics, 
but also undergoes non-unitary evolution as it interacts 
with the environment. For weak system-environment 
couplings, the dynamics of the quantum system is often 
Markovian. However, as the coupling increases, informa- 
tion about the system, which leaks out into the environ- 
ment, may flow back to the system at a later time, mak- 
ing the evolution non-Markovian. Deterministic control 
of the transition from Markovian to non-Markovian dy- 
namics of an open quantum system was recently demon- 
strated in a quantum optical experiment [2j. 

In electronic transport, the theory of open quan- 
tum systems can describe non-equilibrium charge flow 
through nano-electronic conductors. The conductor ex- 
changes particles with the external electronic leads and 
may also interact with a heat bath, Fig.[Il The full count- 
ing statistics of transferred ch arges [3-8] as well as the 
finite-frequency current noise have been investi- 

gated intensively for non-Markovian transport processes, 
and shot noise measurements have revealed strong mem- 
ory effects in the transport through a quantum dot in res 
onancc with the Fermi level of an external electrode 
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A very recent interest in quantum transport concerns 
the distribution of waiting times between consecutive 
charge transfers IB- H| . This line of research seems par- 



FIG. 1. (color online). Generic quantum transport setup. A 
nano-scale electronic system is connected to source and drain 
electrodes as well as an external heat bath. The applied bias 
drives electrons through the system from the source to the 
drain electrode. The system exchanges energy with the bath, 
but no particles are transferred between system and bath. We 
consider situations where the dynamics of the system is non- 
Markovian due to the electronic reservoirs and the heat bath. 



In this Letter we consider quantum transport governed 
by a generic non-Markovian GME and derive a gen- 
eral expression for the electronic WTD in Laplace space. 
The GME contains the memory kernel of the transport 
process and an inhomogeneity term, which accounts for 
memory effects at the time, when charge detection be- 
gins. This term has been the subject of recent theoretical 
investigations flij , 13 1. It typically decays on the same 



ticularly relevant in the light of the increasing number of 
accurate single-electron counting experiments !19| . Theo- 
ries have now been developed to describe electronic wait- 
ing time distributions (WTDs) for driven single-electron 
and phase-coherent conductors 
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quantum optics, non-Markovian effects in the decay dy- 
namics of laser-driven systems have been examined [20] 
and stochastic simulations are currently being used to 
extract the WTDs of non-Markovian quantum systems 
21]. In electronic transport, Brandes has developed a 
compact and elegant method to calculate WTDs for sys- 
tems described by Markovian generalized master equa- 
tions (GMEs) [Hjj. 



time scale as the memory kernel and is not important for 
the long-time limit of time-integrated quantities like the 
zero- frequency noise and higher cumulants of the current, 
but is crucial to include when investigating correlation 
functions and fluctuations at finite times and frequen- 
cies 

Bom- As we demonstrate below, it must also 

be incorporated together with the kernel in a consistent 
theory of WTDs for non-Markovian quantum transport. 

We illustrate our methodology by evaluating the WTD 
for a dissipative double quantum dot (DQD) whose 
dynamics is non-Markovian due to a strongly coupled 
heat bath. This electronic analogue of an open spin- 
boson problem provides us with a microscopic model 
that can describe the transition from Markovian to non- 
Markovian dephasing. We show how coherent oscillations 
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between the quantum dots are washed out by an increas- 
ing bath temperature, which dephases electrons as they 
propagate through the DQD, and we discuss WTDs as 
the system-bath coupling becomes strong. We then tune 
the DQD to a parameter regime, where the WTD be- 
comes particularly sensitive to the spectral properties of 
the heat bath, in a similar spirit to proposals for detecting 
the high-frequency quantum noise of a mesoscopic con- 
ductor by measuring the noise-induced inelastic current 
in a nearby DQD [22| . Within this approach we demon- 
strate how the absorption and emission of energy quanta 
to and from the heat bath can be clearly identified in the 
electronic WTD. 

Non-Markovian GME. — We consider a generic non- 
Markovian GME of the form @, % 

T J(n,t)=T / dt'W(n-n',t-t')p(n',t')+j(n,t), (1) 

describing charge transport through a nano-scale conduc- 
tor as illustrated in Fig. Q] Here p(n, t) is the reduced 
density matrix of the quantum system, obtained by trac- 
ing out the external electronic reservoirs and the heat 
bath. It has been resolved with respect to the number 
of transferred electrons, such that P(n,t) — Tr{p(n,t)} 
is the probability of having collected n electrons in the 
drain during the time span [0,t] (23[. We assume low 
electronic temperatures compared to the applied voltage 
so that thermal charge fluctuations are negligible and we 
may focus on the uni-directional non-equilibrium charge 
transport from source to drain with n being non-negative. 

The kernel W(n,t) determines the time evolution of 
p(n,i), taking into account memory effects due to the 
electronic reservoirs and the external heat bath. The 
inhomogeneity 7(n, t) describes memory effects from be- 
fore t = 0, when charge detection begins. The first term 
on the right hand side of Eq. (JT]) does not include such 
memory effects. Non-Markovian GMEs as Eq. (1) arise 
in a variety of contexts, for example in the real-time dia- 
grammatic techniq ue J24| and in the Nakajima-Zwanzig 
projection method [25]. They have recently been investi- 
gated in connection with full counting statistics and 
finite- frequency noise 0, [H, UH ■ As we go on to show, 
they also provide a useful starting point for calculating 
WTDs in non-Mar kovian quantum transport. 

Distribution of electron waiting times. — To find the 
WTD, we first solve Eq. ([T]) for p(n,t). To this end, we 
employ an operator-valued generalization of the gener- 
ating function technique in Laplace space by introduc- 
ing the transformed reduced density operator p(s, z) — 
YlnLo Jo" dtp(n, t)s n e~ zt and similarly for W(s, z) and 
7(s,z) [26]. Within this framework, the reduced density 
matrix is readily obtained as p(s,z) = Q(s, z){ , y(s, z) + 
p stat }, where Q(s, z) = [z — W(s, z)]~ x is the resolvent of 
the kernel and p stat is the stationary state. In the dis- 
tant past, the system is prepared in an arbitrary state, 



but the stationary state is reached before detection of 
charges begins at t = 0. The stationary state satisfies 
W(l,0)/5 stat = with the normalization Tr{p stat } = 1. 

For stationary processes, the WTD is related to the idle 
time probability II(t) as W(r) = (t)c^II(t), where the 
mean waiting time is determined by the average particle 
current as (r) = 1/(7) [H E^]. Here II(t) is the prob- 
ability of not observing any electrons in a time interval 
of length r, [to, tg + t]. This probability is independent 
of to and we may take to = 0. The average particle cur- 
rent is well-known and reads {I) = Tr{<9 s W(s, 0)/5 stat } s= i 
In Laplace space, the WTD becomes W(z) = 
(t)z[zII(z)-1]+1, where U(t = 0) = 1 and d T U(r = 0) = 
— l/(r) have been used. Importantly, the idle time prob- 
ability can be expressed in terms of the reduced density 
matrix as H(z) = P(n — 0, z) = Tr{p(s = 0, z)}. Insert- 
ing the solution for p(s, z), we arrive at the key result of 
this section 



W(z) = (t)z z(g{z)) -1+1 



(2) 



with (g(z)) = Tr 0(0, z){j{0, z) + p stat } . This equa- 
tion generalizes the expression by Brandes for Markovian 
GMEs [l5| ■ It is easy to show that our result reduces to 
that of Ref. [ill, if Eq. ([1} is replaced by a Markovian 
GME. We see that the WTD corresponding to a generic 
non-Markovian GME includes both the kernel and the 
inhomogeneity 7 To revert the WTD to the time 

domain an inverse Laplace transformation must be per- 
formed analytically or numerically as illustrated below. 

Dissipative DQD. — As a concrete application of our 
method, we consider charge transport through a DQD, 
where non-Markovian effects occur due to a strongly cou- 
pled heat bath. The total Hamiltonian of the setup reads 



H — Hs + Ht + Ht, + Hsb + Hi 



(3) 



where H s = (e/2)(Qd L ~ $ R d R ) + T c (d\d R + aV R d L ) de- 
scribes the left and right levels of the DQD with dealign- 
ment e and tunnel- coupling T c ; tunneling between left 
(right) level and left (right) lead is accounted for by 
H T = J2k,a=L,R t k-acl a d a +h. c; the Hamiltonian of 
the leads is H L = J2k, a =L,R e k,acl a c k , a ; the coupling 
between the DQD and the heat bath reads Hsb — 
(d^dL-d^dR) Sjfe/2)(a] + aj); and H B = J2j f^j^j 
describes the heat bath as an ensemble of harmonic os- 
cillators. This is an open spin-boson problem, where the 
role of the spin is played by the two single-particle lev- 
els of the DQD which are coupled to a bath of bosons. 
It is a transport problem as charges enter and leave the 
pseudo-spin states from the voltage-biased electrodes. 

We derive a non-Markovian GME for the populations 
of the DQ D by tracing out the electronic leads and the 
heat bath [29(. The reduced density matrix p(n,t) — 
[po(n,t), pL(n,t), pn(n,t), pD(n,t)] T contains the proba- 
bilities for the DQD to be empty, having left or right dot 
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FIG. 2. (color online). Electronic waiting time distributions for a dissipative double quantum dot. (a) Coherent oscillations 
between the quantum dots are visible for weak couplings to the heat bath (a = 0.01). The electrons are dephased by an 
increasing bath temperature which washes out the coherent oscillations, (b) Electronic waiting time distributions beyond the 
weak-coupling limit. The increased coupling to the heat bath tends to localize electrons on the quantum dots. 



occupied, or to be doubly occupied. In Laplace space, 
the kernel reads 



(s,z)-- 



( -T L J) sT R \ 

r L -f+(z) r_(z) sT R 
o r+(z) -r_(z)-r L -r fl o 
V o o r L -v R J 



(4) 

where the tunneling rates between the DQD and the leads 
r a (e) = 2irJ2 k \tk,a\ 2 8(e - e fc , Q ) = T a , a = L,R, are 
assumed to be constant. A large bias across the DQD 
ensures that the (broadened) levels of the DQD are well 
inside the bias- window [30( . The inter-dot tunneling rates 

f ± (z) = T 2 [g\{z ± )+g-(z T )] (5) 

are given by the bath-correlation functions g±(t) = 



e -^(Tt) 7 wn ere z± = z ± ie + (T L + T R )/2 and F(t) = 
dLdJ(u){[l-cos{ujt)]coth(huj/2k B T) + ism{ujt)}/uj 2 . 
Here, J(w) = ■ \gj\ 2 6(to— uij) is the spectral function of 
the heat bath with temperature T. We consider an ohmic 
bath with J(uj) = 2auje~ UJ ^ c , where a is the strength 
of the coupling to the DQD and w c is a high-frequency 
cut-off. For strong couplings to the heat bath, the inter- 
dot tunneling rates are valid to lowest order in T c 2 . The 
inhomogeneity reads j(s,z) = — W(l, z)p stat /z 0, El]. 
Strong Coulomb interactions between the quantum dots 
may be included by excluding the double-occupied state 
of the DQD. 



Non-Markovian dephasing. — We evaluate the WTD in 
Laplace space using Eq. ^ and find the analytic result 



W(z) = 



T L T R (z + T L + r R ) 2 r+{z) 



(z + T L )(z + T R )(T L + r^)[zr_(z) + {z + T L + T R }{z + T+(z)}} 



(6) 



In general, the electron waiting time is determined in 
an interplay between the time-scales associated with the 
incoming electrons in the transport window and those 
associated with the nano-structure. The incoming elec- 
trons are on average separated by the mean waiting time 
f = h/eV, where V is the applied voltage [18|. In our ex- 
ample, the voltage is much larger than the energy-scales 
of the DQD, implying that the mean waiting time be- 
tween the incoming electrons is much shorter than the 
time-scales of the DQD, and f does not appear above. 

Focusing first on the uncoupled case (a = 0), we find 
that the WTD in the limit of vanishing tunneling rates 
(J*Li^ R —¥ 0) has imaginary poles at z = ±iA, where 



A = y/ 4T C 2 + e 2 is the energy splitting of the hybridized 



states of the DQD, see also Ref. [15|. These poles corre- 
spond to coherent oscillations between the quantum dot 
levels with period h/A. This is clearly visible in Fig.[2k,, 
showing the WTD in the time domain. The pole struc- 
tures Ti,m/{z + r £ /r) in Eq. (j6|) are due to poissonian 
charge transfers between the DQD and the leads, which 
damp the oscillations. The coherent oscillations are grad- 
ually washed out as the temperature of the heat bath 
is increased and electrons are dephased as they tunnel 
through the DQD. 

Next, we increase the coupling to the heat bath. 
Markovian dephasing in the weak-coupling limit has been 
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FIG. 3. (color online). Electronic waiting time distributions 
for dealigned levels. At low temperatures (blue lines), there 
is a strong asymmetry between positive and negative dealign- 
ments. The asymmetry diminishes at high temperatures (red 
lines). Red and blue curves are obtained from Eq. while 
the black curves follow from the approximation in Eq. ([7|. 



exact results based on Eq. J6|). 

Conclusions. — We have presented a theory of elec- 
tron waiting times for non-Markovian generalized mas- 
ter equations which unifies and generalizes a number of 
earlier approaches to waiting time distributions in the 
context of electronic transport. As an illustrative exam- 
ple we considered electron transport through a double 
quantum dot for which we examined non-Markovian de- 
phasing mechanisms beyond the weak-coupling limit. We 
hope our method may pave the way for future investiga- 
tions of memory effects and electron waiting times, sim- 
ilar to how full counting statistics and finite-frequency 
noise in non-Markovian quantum transport have been 
popular research topics in recent years. 
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discussed by Brandes [15|. Here we can take these ideas 
further and examine the transition from Markovian to 
non-Markovian dephasing. In Fig. [3Jd, we present WTDs 
at zero temperature beyond the limit of weak system- 
bath coupling. As the coupling increases, the heat bath 
tends to localize electrons on the quantum dots and the 
inter-dot tunneling rate becomes suppressed. For large 
couplings, tunneling events are rare and uncorrelated and 
the transport process essentially becomes poissonian. 

Heat bath & WTD. — The DQD can be tuned to an 
interesting regime, where tunneling between the quan- 
tum dots becomes the rate-limiting step in the trans- 
port. Choosing for instance the tunnel coupling T c or 
the dealignment of the quantum dot levels e such that 
r_(z),r + (z) < T L ,T R , Equation © reduces to 



W(z) 



z + r+(z)' 



(7) 



This result offers the possibility of directly probing spec- 
tral properties of the heat bath (or another coupled con- 
ductor [32|) through the detection of the electron wait- 
ing time, since the bath correlation functions enter the 
bath-assisted hopping rate T + (z). In Fig. [3] we focus 
on the emission and absorption of energy to and from 
the heat bath as electrons tunnel from the left to the 
right quantum dot. At low temperatures, there is a clear 
asymmetry between the WTDs for positive and negative 
dealignments, since the heat bath mainly contributes to 
the transport for positive detunings by absorbing energy 
quanta from tunneling electrons. At high temperatures, 
this asymmetry disappears as the heat bath in addition 
can assist the tunneling process at negative detunings 
through the emission of energy quanta. Figure 3 shows 
that Eq. ([7]) provides an excellent approximation to the 
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